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Abstract: In direct dark matter detection experiments, conventional elastic scattering 
of WIMPs results in exponentially falling recoil spectra. In contrast, theories of WIMPs 
with excited states can lead to nuclear recoil spectra that peak at finite recoil energies E^. 
The peaks of such signals are typically fairly broad, with AEn/Epeak ~ 1- We show that 
in the presence of dark matter structures with low velocity dispersion, such as streams 
or clumps, peaks from up-scattering can become extremely narrow with FWHM of a few 
keV only. This differs dramatically from the conventionally expected WIMP spectrum and 
would, once detected, open the possibility to measure the dark matter velocity structure 
with a fantastic accuracy. As an intriguing example, we confront the observed cluster of 
3 events near 42keV from the CRESST commissioning run with this scenario, and find 
a wide range of parameters capable for producing such a peak. We compare the possible 
signals at other experiments, and find that such a particle could also give rise to the signal 
at DAMA, although not from the same stream. Over some range of parameters a signal 
would be visible at xenon experiments. We show that such dark matter peaks are a very 
clear signal, and can be easily disentangled from potential backgrounds, both terrestrial or 
due to WIMP down-scattering, by an enhanced annual modulation signature in both the 
amplitude of the signal and its shape. 
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1. Introduction 

The search for dark matter is inextricably bound up with our theories of what it might 
be. With the only robust signals for its existence being gravitational in nature, we are left 
to speculate from outside motivations as to its nature and hence on how we might find 
it. Since search strategies depend intimately on precisely those ideas, it requires a delicate 
blend of astrophysics and particle physics, with assumptions from both, to produce a 
final signal for which experiments can search. To this end, particle physics uncertainties 
can be studied through the development of new and varied models, whereas astrophysical 
uncertainties can be studied through variations in halo models, and more recently through 
the application of realistic N-body simulations. 

In the context of conventional WIMPs, signals from astrophysics are somewhat limited 
for most direct detection experiments. While streams can change the phase and amplitude 
of a modulated signal [1], and subhalos can dominate the signal at a directional experi- 
ment [2], for the majority of currently and soon-to-be operating experiments, the signal is 
not profoundly different. For spin-independent interactions in particular, the exponentially 
falling nuclear form factor compresses events to the lowest energies, yielding quantitatively 
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different - but qualitatively similar - signals, at least within the next round of experi- 
ments [3]. 

At the same time, a wide variety of models have arisen of late that are sensitive to 
the motions of particles at the highest velocities in the halo. Inelastic dark matter [4, 5], 
resonant dark matter [6], and dark matter with momentum-dependent scattering [7, 8] all 
can yield signals at high energies where the sensitivities to structures of the halo become 
important. A recent study [2] employing the Via Lactea II and GHALO simulations has 
shown that, in the particular case of inelastic dark matter, the variations of the halo can 
change the spectrum and sensitivities of the various experiments significantly. 

In the inelastic dark matter scenario, dark matter scatters off of nuclei by transitioning 
to an excited state an amount S ~ m^fi'^ heavier than the ground state. These kinematical 
changes can lead to dramatic changes in the expected signal: a change in the relative 
sensitivity of different nuclei, an enhancement of the annual modulation and, perhaps most 
importantly from the perspective of a direct detection experiment, a peak of the signal at 
relatively high energies, with few or no events at lower energies. 

In the context of the recent PAMELA measurement of an excess of positrons at high 
energies [9], such new signals receive new and independent motivation. To explain the 
large annihilation cross section required by PAMELA with thermal WIMPs, one must 
invoke a new force carrier to yield a low velocity increase through the Sommerfeld en- 
hancement [10]^ or capture into WIMPonium [15]. The light force carrier is required to 
produce a large enhancement, but additionally naturally produces copious hard leptons 
without any antiprotons [16] consistent with observations [17, 18]. If this new mediator is 
a gauge boson, the WIMP must automatically contain multiple states, either multiple real 
scalars or multiple Majorana fermions. For the Sommerfeld enhancement to be effective, 
those states cannot be significantly heavier than a^m^^, or else the annihilation rates would 
be suppressed [19]. This suggests that not only should the excited states be there, they 
would also likely he kinematically accessible, modifying our expectations for direct detection 
experiments, irrespective of the outside motivation from DAMA. 

In this paper we explore in more detail the effects of the substructures in dark matter 
halos for direct detection experiments. We begin by comparing the elastic and inelastic 
cases in section 2, and detail how inelastic scattering in the presence of cold streams 
or subhalos can lead to sharply peaked signals. In section 3, we consider in particular 
the CRESST experiment, which saw seven events in its commissioning run. We describe 
what ranges of signals would fit the events observed in section 4 and argue that seasonal 
effects should be 0(1) for substructure signals, allowing one to positively identify them. 
Specifically, this will allow a discrimination between a heavy WIMP upscattering and a 
light WIMP down-scattering (section 4.4). In section 5, we consider the implication of 
these scenarios for other experiments such as CDMS or XENON and study the effects of 
streams on DAMA. 



The importance of this effect in the context of multi-TeV scale dark matter interacting via W and Z 
bosons was first discussed by [11], and more recently emphasized with regards to PAMELA in the context 
of "minimal dark matter," with similar masses, by [12, 13, 14]. 
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2. Recoil Spectra from Dark Matter Velocity Substructure 



2.1 Elastic Scattering 

The standard signal expected at a direct detection experiment is tiiat of a WIMP elastically 
scattering off of a nucleus, yielding an exponentially falling distribution (figure 1). Let us 
take a moment to understand this spectrum from spin-independent scattering, which will 
help us understand the sensitivity of the spectrum to velocity substructure in the case of 
inelastic scattering as well. 
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Figure 1: The spectrum of an elastically scattering WIMP. Dashed (red): A monoenergetic stream 
of elastically scattering WIMPs produces a flat spectrum up to a maximum energy fi'max- Solid 
(blue) : taking coherence loss into account by means of the form factor strongly suppresses the high 
energy recoils, in particular for heavy target nuclei. Dotted (grey): Taking the velocity distribution 
of WIMPs into account (here assumed to be of Maxwcll-Boltzmann type) gives the usual quasi- 
exponential spectrum. 



The standard nuclear recoil spectrum for spin-independent scattering is given by 



^target "I AT 



^-2<^A'F\En) r dv^, (2.1) 



where A^target, ttT'N and A is the number, mass and mass number of the target atoms, 
and the WIMP mass and density, an the zero-momentum spin-independent WIMP- 
nucleon scattering cross section, /ine is the reduced mass of the WIMP-nucZeon system, 
the nuclear form factor, fmin the minimum WIMP velocity needed to create a recoil with 
recoil energy Er, and / the WIMP velocity distribution. We assume the couplings to all 
nucleons are equal, for simplicity. 

Let us consider the scattering of a population of WIMPs moving at a single velocity 
■Wstrcami i-S- f{v) = 5{v — ^stream)- In this case, the above equation simplifies to 
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Hence, the expected spectrum is flat before taking into account the nuclear form factor 
effects (see figure 1). The upper hmit of this distribution is simply set by the maximum 
energy transfer from the WIMP to the nucleus, 



where rriAr is the mass of the target nucleus and Vj, = t'wiMP/^/wiAr, with fi here the 
reduced mass of the WlMP-nucleus system. The nuclear form factor F'^{Eji) is in general 
an exponentially falling function, which suppresses higher energy recoils, yielding a falling 
spectrum. Furthermore, after one integrates over a reasonable velocity distribution f{v), 
with fewer high energy particles capable of scattering to the highest Eji, this is accentuated. 
The combination is the well-known nuclear recoil distribution of WIMP elastic scatters, 
which falls monotonically until the first form factor zero, after which it is generally very 
small. 

Extracting information about the halo from these exponentially falling spectra is very 
challenging [20, 21], and for elastically scattering WIMPs, detecting velocity substructure 
seems to be nearly impossible. In general, standard WIMP induced spectra are not par- 
ticularly sensitive to velocity substructure since they are caused by WIMPs with a wide 
range of velocities. Hence, any contribution from substructure is dwarfed by the smooth 
halo component. 

To detect substructure, at least straightforwardly, one must appeal to models in which 
a narrower range of WIMP velocities can contribute to the spectrum at a given energy. 
As we shall see, in the context of inelastic models, sharply peaked signals can arise at 
experiments with good energy resolution, and these can inform us about WIMP velocity 
structures. 

2.2 Inelastic Dark Matter 

A simple scenario which is extremely sensitive only to the high energy tail of the WIMP 
velocity distribution is the inelastic dark matter scenario [4, 5]. In this scenario, it is 
assumed that dark matter x has an excited state x*^ split by (5 ~ where is the 
reduced mass of the WIMP-nucleus system. Let us suppose that the WIMP can only 
scatter by transitioning from the ground state to an excited state, i.e., x^ ~^ X*^: while 
elastic scattering is suppressed. Models for inelastic dark matter are simple to construct [4, 
5, 10, 22, 23, 24, 25, 26], and can have implications for astrophysical signals from their 
capture in compact objects [27, 28, 29, 30, 31]. We shall not address these issues here, 
limiting ourselves to the direct detection possibilities. 

The assumption that the splitting is comparable to the available kinetic energy is 
essential to this scenario, because this significantly alters the scattering process. The 
principle kinematical change is that the incident velocity of the WIMP must be large 
enough to cause excitation. Specifically, the WIMP velocity must exceed 
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It is clear that a high momentum transfer is favored for heavy target atoms, and indeed, 
function 2.4 has a minimum at 
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(2.5) 



The simple change in this scenario is enough to allow a positive signal at DAMA while 
being consistent with other experiments [32, 33, 26]. Let us examine in detail the spectrum 
of such a scenario. 

If we ignore for a moment the nuclear form factor and consider the case of WIMPs with 
a single velocity /3min that is just enough to excite the WIMP, scattering is at threshold. In 
the center-of-mass frame, the WIMP has just enough energy to interact with the nucleus, 
so after scattering, both will cling together. Boosting this scenario in the lab frame gives 
the nucleus the velocity of the two frames V, which is given simply by the velocity of 
Earth in the galactic halo and the reduced mass. Thus, there is a non-zero signal only for 
Er = En, so the signal is completely monochromatic. Now, suppose the WIMPs have a 
higher velocity. In the center-of-mass frame, this will leave the nucleus with some velocity 
V, which needs to be added vectorially to V to boost in the lab frame. The recoil energy 
Er can thus be either higher or lower than Er depending on the angle between v and 
V. Therefore, the upper bound of the spectrum is increased and the lower bound of the 
signal is decreased (see figure 2). Before the inclusion of the form factor and any velocity 
dispersion, the spectrum is thus centered around Er and flat. After applying the form 
factor, this flat distribution becomes sharply peaked. 
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Figure 2: The spectrum of an inelastically scattering WIMP. Dashed (red): For a monoenergetic 
WIMP stream, no recoils take place below a given threshold since energy is required to excite 
the WIMP. Solid (blue): The form factor suppresses the higher energy recoils, yielding a strongly 
peaked recoil spectrum. Dotted (grey): Expanding the velocity distribution of the WIMPs to a full 
Maxwell-Boltzmann distribution will smear the peaked recoil spectrum. 



This point is one of the central of our paper: the presence of a population of inelastically 
scattering WIMPs with a low velocity dispersion (i.e., a stream) can lead to sharply peaked 
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signals, with a peak at the energy which is the lowest energy kinematically allowed for 
those WIMPs. Although the extreme example of a stream with zero velocity dispersion 
seems somewhat artificial, adding a small velocity dispersion will simply broaden the peak. 
Moreover, the presence of such narrow peaks has already been found in the application of 
N-body simulations to the inelastic dark matter scenario [2]. 

Before moving to a discussion of these peaks in existing experiments, it is important 
to note that there is a large degeneracy between the parameters 5, my_ and f stream- The 
central value for the distribution (before applying the form factor) is at = Sfi/niN, 
with a width of 



Thus, the observation of a peak will only constrain one combination of these, and the 
measurement of the high and low part of the spectrum will constrain two. In the limit 
of /u — )■ niN, these two will determine ^stream and 6, as drops out of the equations. 
However, one cannot tell if one is in this limit, and thus multiple experiments will be 
required to determine these parameters, as we describe in section 5. 

Also, it is all possible that inelastic dark matter interacts with standard model matter 
in a spin-dependent fashion. Such a possibility can alleviate constraints from direct detec- 
tion experiments [34], in particular CRESST, of whose tungsten isotopes only ^^^W has an 
unpaired neutron. Peaked signals can arise in this scenario, although less so, because the 
spin-dependent form factor does not in general fall as rapidly as the spin-independent. We 
will focus on the spin-independent case, but it is important to note that these signals are 
more general. 



The CRESST-II experiment [35, 36] fulfills the requirements to measure these parameters 
in an excellent way. It operates CaW04 crystals at temperatures of ~ 15 mK, with tungsten 
being the heaviest target used in any direct detection experiment today. In addition, the 
experiment achieves an excellent energy resolution of AE < 0.6 keV (Ifi) for energies E 
below 100 keV, reaching AE = 130 eV (la) at E = 3.61 keV [37]. 

The energy scale quoted by the collaboration is in terms of keV electron equivalent 
(keVee), inferred from a calibration with 122 keV gammas from ^^Co. Quenching effects 
are neglected, which can be justified in the search for the usual WIMP signatures, since 
CRESST detectors are good calorimeters [37] and at most 5% of the total interaction 
energy is converted into scintillation light [38] . Here, however, we are dealing with strongly 
peaked signatures and need to take quenching into account. From [39], one can estimate 
that 1.3% of the interaction energy are detected as scintillation light, but needs to factor 
in light that escapes the crystal but is not detected. Hence, we estimate that the energy 
measured in terms of nuclear recoil equivalent keVR is 3% lower than the value in terms of 
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In fact, the experiment has observed nuclear recoils that resemble the peaked signature 
expected from inelastic dark matter streams. We consider data taken after a neutron 
shield was installed to reduce the neutron induced nuclear recoil background [40], but 
stress that data taken earlier [35, 41] are consistent with our observation. A peak at 
42keVee ~ 41keVR stems from three events that fall within 0.7 keV of each other. The 
chance probability for such a cluster of events can readily be estimated to be a few percent 
at most, given the low apparent event rate. However, the distinct signature discussed here is 
difficult to create otherwise. The dominant neutron induced nuclear recoil background is a 
simple falling exponential just like the usual WIMP spectrum [42], but otherwise featureless 
and hence very different from the peaked dark matter spectrum discussed here. In addition, 
in the CaW04 target of the CRESST-II experiment, neutrons dominantly cause oxygen 
recoils which can be distinguished from WIMP induced tungsten recoils based on their 
light yield [36]. 

Another possible nuclear recoil background comes from surface alpha decays. The 
most well studied example is that of ^^^Po (a ^^^Rn daughter), which alpha decays into 
206p]-, a half-life of 138 days and a total energy of 5407 keV available in the decay. The 
ejection of the 5304 keV alpha particle is of no concern here, but the lead nucleus recoils 
with an energy of 103 keV and can mimic a dark matter induced nuclear recoil signal [41]. 
If the ^^''Po has been implanted in a surface close to the target by its mother decays, the 
recoiling ^'^^Pb looses energy as it escapes this surface. This can result in a recoil signal at 
lower energies than 103 keV [43]. The spectral shape of these ^^'^Po events is thus peaked 
at 103 keV with a tail extending to lower energies. However, no alpha decay exists that 
would lead to nuclear recoils with energies in the range discussed here. 

While this could be simply a few percent chance alignment, it serves as a useful example 
that we will pursue for the remainder of the paper, as well as its associated questions. What 
ranges of parameters could generate a peaked signal within the experimental ranges of 
existing experiments? What are the sensitivities of other experiments to it? If it persisted, 
how could we determine its dark matter origin? Even if this coincidence proves ephemeral, 
it will provide a useful case study to address these questions. 

4. Peaked Signals from Substructure 
4.1 Generic Considerations 

It is quite simple to find a set of parameters that would work for an inelastic WIMP to 
explain the observed events in CRESST. We assume that the additional subhalo is moving 
directly counter to the sun's motion in the galaxy. To the extent that such subhalos are 
detectable, this is likely not a bad approximation. Halos from directions far from this will 
not gain the boost from the Earth's average motion, and so will generally be dwarfed by 
the smooth halo. This is borne out by the study [2] which showed that the hottest spot 
on the sky (i.e., the region of the sky with the largest flux of WIMPs) for particles with 
■Wstream > 500km/s is generally within 30° of the Earth's motion. It is impossible to state 
the likelihood of such a structure to exist in the local neighborhood. Numerical simulations 
have yet to resolve sub-parsec scales, but because the scenarios we consider are sensitive 
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to the highest velocities only, the likelihood that such a structure could be relevant is 
significantly enhanced compared to a standard WIMP scenario, where a broad range of 
velocities can contribute. Even if the total mass in the stream or subhalo is sub-percent of 
the local density, it can still be the dominant contributor to the signal, because of its high 
velocity, where the smooth component is exponentially suppressed. We shall not further 
address the likelihood of this, but will proceed in attempting to understand what such a 
structure would look like, should it appear. 

In figure 3, we show a simple set of parameters that give rise to a narrow spike at 
~ 41 keV. In particular, a 100 GeV WIMP with 5=130 keV can give such a peak in the 
presence of a subhalo moving towards us with velocity 370 km/s in the halo rest frame, 
corresponding to about 600 km/s in the Earth frame. Also, a 700 GeV WIMP with 
6 = 150 keV can give a similar peak. In this latter case, additional signal is expected at 
higher energies, although the amplitude of this higher energy signal can be suppressed if 
the mediator is < 100 MeV [44]. In figures 3-9 we take = 100 MeV for the spectra of the 
heavier WIMP to suppress this signal somewhat, and m| ^ for the lighter WIMP, where 
it is not present. Searches and low-energy constraints for such particles are summarized 
in [45]. 
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Figure 3: Peaked signals arising from WIMPs (red) and events observed by CRESST (blue) in 
the commissioning run [40]. Left: = 100 GeV, S = 130 keV, and a cold stream with velocity 
370 km/s in the galactic rest frame (^ 600 km/s is the lab frame). Right: rii^ = 700 GcV, S = 
150 keV, and stream velocity 280 km/s (~ 510 km/s). 



Of course, CRESST observed more events than just at 41 keV, and while they may be 
background, it is worth considering whether they could be arising from the same inelastic 
WIMP. Not surprisingly, it is quite easy to fit them. In figures 4 and 5, we consider the 
same model as in figure 3, but with an additional Maxwellian component or additional 
subhalos, respectively. 

4.2 Seasonal Variation 

The question arises how such peaked, inelastic dark matter induced signals could be dis- 
entangled from terrestrial backgrounds, given the available experiments. If the experiment 
that makes the positive detection uses a fairly light nucleus, such as germanium or lighter. 
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Figure 4: As in figure 3, but with an additional MaxwcU-Boltzmann contribution in addition to 
the stream. 
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Figure 5: As in figure 3, but with additional streams. Left: m-^ = 100 GeV, S = 130 keV, and a 
cold streams with velocities 370, 385, 450km/s in the galactic rest frame. Right: m-^ = 700 GeV, 
S = 150 keV, and stream velocities 280, 310, 425km/s. 



then the xenon experiments (XENONIOO, XMASS, LUX) will be kinematically capable 
of detecting the same structure and the multi-target approach allows to disentangle signal 
from background. On the other hand, if the experiment with a signal uses a heavy nucleus, 
such as tungsten, then it may well be that no other experiments are kinematically capable 
of detecting the same WIMP. In that case, testing for an annual modulation of the sig- 
nal within the same experiment will give a clear indication. The annual modulation can 
give an indication of substructure, while its phase and amplitude can give evidence for the 
direction of the streams motion. 

Let us begin by considering the case considered so far - that of a stream which is 
moving directly counter to the galactic rotation. Due to the significant inclination of the 
Earth's orbit relative to the galactic plane, this is not the situation for maximal modulation. 
We assume that the signal is arising from a cold stream near the galactic escape velocity. 
Although the change of velocity distribution from summer to winter is only 0(30 km/s), 
this is extremely significant. We show in figure 6 the summer and winter spectra for a 
heavy and light WIMP case where a peak contributes the only visible signal, here at 41 
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keV. We take the summer to be ±45 days around the peak, and winter ±45 days around 
the minimum. It can be seen that for a hght WIMP, the entire signal disappears in winter, 
while for a heavy WIMP, the peak disappears and the signal at high energies persists. 




£„/keV E„/keV 

Figure 6: Annual modulation of the signals considered in figure 3. Solid (blue): Signal in summer, 
dashed (red): Signal in winter. Left: = 100 GeV, S = 130 keV, ^stream = 370km/s in the 
galactic frame. Right: niy, = 700 GcV, 6 = 150 kcV, ^stream = 280km/s . 



The modulation for these peaks can be very brisk. For the case of the 41 keV peak, 
in particular, the signal can vanish rapidly. We show in figure 7 the amplitude in the peak 
as a function of time. In the case of the lighter WIMP, the signal is only present in May 
and June, while in the heavier case, it is present for a longer period, but is suppressed 
significantly still by the end of August. 
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Figure 7: The signal present in the range 39-45 kcVr as a function of month (0=June 2). Solid: 
= 100 GeV, 5 = 100 keV. Dashed: = 700 GeV, 5 = 150 keV. 

One can compare this with he situation arising from a Maxwellian distribution, which 
we show in figure 8. In contrast to the rapidly varying peak, the Maxwellian structure 
largely retains its shape and has 0(20%) modulation. 

In the presence of multiple peaks, generally the highest energy peaks will be diminished 
or lost, while the lower energy peaks will remain, which we show in figure 9. As seen 
before, the high energy signal above the form-factor zero is generally constant throughout 
the year. Below the form-factor zero, however, spikes can shrink dramatically or vanish 
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Figure 8: As in figure 6, but with an additional Maxwell-Boltzmann component, with relative 
normalization as in figure 4. 



entirely, depending on the parameters. Thus, even with only a single experiment, the 
dramatic annual modulation should arise with adequate statistics, and give confirmation 
that a dark matter substructure signal has been observed. 
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Figure 9: As in figure 6, but with additional streams. Left: Streams at 385 and 450km/s in the 
galactic rest frame. Right: Streams at 310 and 425km/s, with relative normalizations as in figure 
5. 



4.3 Amplitudes, Phases and the Sagittarius Stream 

Up to this point we have limited ourselves to the case of a stream moving in the motion 
directly opposite that of the Galactic rotation. This leads to the highest possible velocities 
in the local frame, but as the velocities we have considered have been well below the 
Galactic escape velocity, one can achieve similar effects with streams moving in other 
directions. In particular, the most relevant question is the direction of the stream in the 
solar frame, relative to the Earth's motion. This will determine the amplitude and phase 
of the modulation, which are detectable parameters. Should the motion of the stream be 
directly opposing the Earth's motion at some point in the year in the solar frame, the 
modulation can be even larger than we have discussed. If the direction is not opposite the 
Earth's motion in the summer, the phase can shift dramatically. 

One possible source of a cold substructure in the solar neighborhood is the Sagittarius 
stream, see e.g [46]. The possible implications in the context of dark matter searches have 
been well explored [47, 48, 49, 1], in particular in the context of light WIMPs. It is not 
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clear whether the stream is even present in the local neighborhood, see the recent discussion 
in [50]. However, it serves as an example, and even if stellar debris seems absent locally, 
an extended dark matter structure might still be present. What sort of signal might it 
produce? 

The Sagittarius stream has been studied and modeled extensively [51, 52, 53, 54, 46, 55]. 
To model the Sagittarius stream, we employ the results of the simulations by [46], using in 
particular the results of the oblate simulation, where the stream had the highest velocity 
and thus is most relevant for our study. We take particles from the simulation within 
2 kpc of the solar location and find a mean velocity (in the galactic frame) of 415 km/s, the 
mean components of the velocity vector as (72, 75, —400) km/s with norm 413 km/s and 
dispersion 12 km/s, which is an adequately narrow range to treat as a single velocity for 
our purposes. 

Since the Sagittarius stream falls onto the galactic plane roughly from above, its signal 
will be strongest in winter, with peaks appearing in winter but vanishing in summer. It 
seems difficult at best to explain the 41 keV cluster at CRESST, as those data were taken 
in the summer [40] , and previous runs in the early spring (which had no neutron shield and 
thus higher backgrounds) did not show any dramatic excess [35]. Still, the appearance of 
a spike peaking in the winter, and consistent with f stream ~ 400 km/s could be indicative 
of this structure. 

4.4 Peaked Signals from Down-Scattering 

Another possibility for a peaked structure that might fake substructure is inelastic down- 
scattering, for instance from metastable WIMP states [56, 57]. Let us consider a WIMP 
incident upon a nucleus with velocity v, which scatters from an excited state x* down to 
the ground state x- Then the energy En of the recoiling nucleus with mass m^r is bounded 
by the incident WIMP kinetic energy 

1 2 1 2 

-niN {vb - Vr) <Er< -rUN {Vb + Vr) (4.1) 

where Vb = fi/mi\iv, and Vr = Vb\/l + 26/fiv'^. This leads to an overall width for the line 
of 

2mNVbVr = 2fi^v'^/mN^/l + 26/fiv'^. (4.2) 

If 5 ^ IJ,v'^, this would be indistinguishable from elastic scattering. We are interested in 
the case where 6 ^ fiv'^ and the line would be narrow. In this case, the signal is peaked 
with energy Er = (fiv'^ + 6)n/m]^ 6fi/mN and width 6E/ Er A\/ ^v"^ /25. Hence, the 
line is narrow when the characteristic energy of the center-of-mass system is smaller than 
5. If we are interested in peaked signals in the range Er < 100 keV, we are compelled to 
consider two options: either the velocity of the WIMP is small, or the mass of the WIMP 
is small. 

Although there is evidence from some numerical simulations that there could be a 
dark disk with a low velocity dispersion [58], the velocity dispersion is still expected to be 
~ 30 km/s. Even if it were smaller, it would need to be co-rotating with the sun, whose 
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motion relative to the local standard of rest is significant [59] and includes a sizable motion 
~ 7km/s perpendicular to the plane of rotation. Thus, a heavy WIMP with low velocity 
down-scattering will likely not give a sharp feature. On the other hand, assuming a typical 
velocity of ~ 220km/s, one can achieve a ~ 10% width for down-scattering WIMPs with 
masses < 1 GeV. 

Can we distinguish such a scenario from inelastic substructure? In fact, it should be 
quite straightforward. First off, the down-scattering signal has very little modulation, in 
particular in its shape. Secondly, a down-scattering signal should appear on every target, 
subject to an adequate exposure. Moreover, the peak should shift in a well defined fashion 
for the different targets. In the limit that we take <C mi\f for all experiments, the 
location of the peak should be at -Epcak oc l/rriN- Thus, a simple comparison between 
signals at different targets should clarify this origin. 

As an example, we can ask how a down-scattering signal would appear for the 41 keV 
CRESST cluster. Taking a 1 GeV WIMP, fixing Sji/mM = 41keV, and normalizing to 
three events at CRESST, we would expect > 20 events on germanium targets with current 
exposure (e.g. CDMS and EDELWEISS), and ~ 8 events at XENONIO. For XENONIO, 
the events would lie at AlMy/ /Mxe = 56keV. A quick glance at [60] suggests that this 
is not the case, with only two events in the vicinity. However, if the employed scintillation 
efficiency Leg = 0.19 is suitably far off, then it is possible the events at lower energy would 
be reinterpreted, and so consistency might be possible. The events at CDMS would be 
peaked at ^ AlMw/Mce = 101 keV, which is outside the WIMP search range of CDMS. 
EDELWEISS has reported data on 144 kg day exposure up to 200 keV, with no events at 
high energies [61]. With an expected signal of more than 20 events (even making conserva- 
tive assumptions about the tungsten form factor) this is clearly excluded. Should a future 
peak arise, such a signal would not modulate significantly, and so a sizable modulation in 
the CRESST data would indicate a more traditional inelastic dark matter scenario. 

5. Signals at Other Targets 
5.1 Implications for DAMA 

The DAMA annual modulation signature has persisted and grown to 9.3(T [62]. Inelastic 
dark matter was proposed as an explanation for this signal in the context of constraints 
from CDMS and later XENON, ZEPLIN, KIMS and CRESST. Comparing the different 
experiments is challenging, because they are generally sensitive to different components 
of the halo, and normalizing to DAMA requires a subsequent extrapolation to different 
components of the velocity distribution. This is important, because the non-Maxwellian 
profiles can alter these relative components significantly [63, 33], and substructures in 
particular can change the amplitude and spectrum of signals [2] . 

Here we can address a limiting case by considering signals only due to a stream. This 
allows a sort of apples-to-apples comparison, because it addresses the real question: if 
a given set of particles are visible at DAMA in a given model, what is the implication 
of those same particles at other experiments? In a sense, this provides us with a set of 
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constraints on the inelastic dark matter explanation of DAMA that are nearly independent 
of astrophysical models. 

We show in figure 10 the signal and constraints from various experiments. For the light 
WIMP case (m^ = 100 GeV), we consider a stream moving with velocity 525km/s in the 
galactic frame, with a motion such that it is directly opposite the Earth's summer motion in 
the lab frame. This direction gives the maximum modulation, with other streams giving less 
(although still large) . For the heavy WIMP case (m-)^ = 1 TeV) , we consider a stream with 
velocity 450 km/s in the galactic frame, with a direction that is opposite the Earth's summer 
motion in the lab frame. As in [2], we sum allowed regions for iodine quenching factors from 
0.06 — 0.09 and use the datasets employed [64, 65, 66] as described there, with the most 
recent CDMS results added [67]. Unlike [2] we employ the maximum gap method for all 
experiments, and limit our analysis to the published CRESST commissioning run data [40]. 
It was argued in [2] that one must be wary of using Yellin's analysis techniques [68] in 
situations where the energy resolution is good, simply because substructure can alter the 
spectrum considerably. But in this case, the spectrum is well known, since, by assumption, 
it arises from a single stream of particles. 
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Figure 10: Sample of the allowed ranges of parameter space to explain DAMA in the presence 
of a stream. Left: Wstroam = 525 km/s, Right: ^stream = 450 km/s. Curves are as follows: Thin, 
short-dashed: CRESST commissioning run; thin, long-dashed: ZEPLIN-III; thick, long-dashed: 
ZEPLIN-II; thick, short-dashed: XENONIO; thick, sohd: CDMS-Ge. 



As seen in figure 10, in this "model independent" approach, a signal consistent with 
DAMA is allowed by the other experiments for both small (100 GeV) and large (ITeV) 
masses. The "light inelastic window" which is present with ion channeling [69] is still 
allowed [34] irrespective of these arguments. Taking p stream / Po = 10"^ or 10^^, we require 
a cross section in the range of lO-^^^cm^ < cr„ < lO^^'^ cm . The lower range could arise 
from Z-exchange [4, 5, 23], while the entire range could naturally arise from the exchange 
of a light mediator [10, 70, 71]. 

How does the DAMA signal compare to the CRESST spectrum in this case? Taking the 
same stream and parameters needed to explain the DAMA signal from figure 10, we show 
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in figure 11 the predicted spectra for a low and high mass WIMP case at CRESST (using 
the commissioning period), as well as the same parameters with a Maxwell-Boltzmann 
halo, normalized to the same overall rate. We can see why the limits are now weaker from 
CRESST than with a Maxwell-Boltzmann halo. First, the large modulation associated 
with a stream yields a lower overall rate at CRESST. Secondly, the expected spectra at 
CRESST can be in a much narrower range (which is roughly coincident with events that 
are observed). In essence, a Maxwellian halo generates signal at CRESST from particles 
that do not give rise to signal at DAMA, and at a broader range of energies. Yellin- 
type analyses using Maxwell-Boltzmann halos can then overstate constraints. In contrast, 
the application of our stream technique gives a limit nearly independent of astrophysical 
uncertainties. Specifically, we expect a peak at CRESST in the summer in the energy 
range 15 ~ 30keV. Higher peaks would be possible if the DAMA quenching factor is much 
smaller than expected. Additional signals can arise from other halo components that do 
not contribute to DAMA, naturally. But, for the moment at least, exclusions from other 
experiments (e.g., [72, 34]) are highly halo- model dependent. 
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Figure 11: Spectra at CRESST from a stream that could generate the DAMA amiual modulation 
signal, in counts/keV/kg/day (dru). Thick (blue): ~ 100 GcV, 5 = 100 kcV, solid - stream with 
velocity 525km/s in the galactic rest frame; dashed - Maxwellian velocity distribution. Thin (red): 
= 1 TeV, S = 135 keV, solid - stream with velocity 450km/s; dashed - Maxwellian distribution. 




5.2 Xenon and Germanium Targets 

It is possible for CRESST to see signals that other available experiments can not, since 
CRESST employs the heaviest target material and has a low energy threshold. However, it 
is worth considering whether any of the scenarios described above would lead to dramatic 
signals on other targets, and where those might arise. We can group available experiments 
with sufficient sensitivity into intermediate mass targets (xenon and iodine) and light tar- 
gets (germanium). Targets much lighter than germanium will find difficulty detecting any 
signal given the splittings we have considered here. 
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Let us focus on the possible signals at a xenon experiment. For a given {m^,5}, 
there is a specific ^stream which will generate a spike at a specific energy Er, which we 
take to be ~ 41 keV. One can ask whether that same set of parameters will generate 
any signal on a different target (i.e., whether ^stream > i'min(Xe)). We illustrate this in 
figure 12 and note that much of the existing parameter space is visible only at CRESST. 
However, there are regions both at low and high 6 also visible at xenon. These regions 
are characterized by peaks as well, although smeared by the energy resolution of these 
experiments. Nonetheless, for instance at high m^^. and high 5, a peaked signal would also 
appear at a xenon experiment. Thus observing both peaks would constrain us to lie on a 
line in — 6 space. 

This is achieved simply by observing the location of the low-energy peak, but higher 
or higher 5 parameters that achieve such a peak tend to also have a large tail extending into 
higher energy (see figure 3). While the overall amplitude of this is subject to uncertainties 
in the nuclear form factors at large momentum transfer, as well as effects from a light 
mediator (p in the propagator (which suppresses the signal if ~ = 2mNEji), the high 
energy cutoff also places kinematical constraints on the parameters. The observation of 
this high energy tail at one (or more) experiment can then completely constrain (and even 
over-constrain) the parameters, determining '^^streami 

and 6. 

While we have focused on xenon, one can repeat this exercise as well for a germanium 
target. For the parameters we have chosen here (namely, a 41keV peak on a tungsten 
target), only a small region of the parameter space at low rny.,6 would be visible on ger- 
manium. However, a peak at lower energies on tungsten (and thus arising from a higher 
velocity stream) could be visible on future germanium experiments as well. 

Moreover, this discussion applies only to the peak at 41 keV and the velocity structure 
that generates it. If this structure is accompanied by a smooth halo component (or an 
additional WIMP state which is closer to the ground state in mass), then signals would be 
expected at lighter targets as well. Indeed, with a galactic escape velocity of ~ 600km/s, 
a signal on a Xe target should appear for essentially the entire parameter space, although 
with unknown normalization. Moreover, a comparison of the parameters in figure 12 with 
those in [2] shows that much of the interesting range is compatible with what is required 
to explain the DAMA annual modulation signal. While the results are logically distinct, 
should the 41 keV spike persist in future CRESST data, it would give support to the 
inelastic dark matter interpretation of DAMA. 

6. Conclusions 

Signals of elastic WIMPs at direct detection experiments generally fall monotonically with 
energy, making an extraction of WIMP and/or halo parameters a challenge. For inelastic 
dark matter in contrast, peaked signals are expected. 

We have explored in detail how such peaks arise, and, in particular, how very narrow 
peaks can arise in the presence of subhalos or cold streams. We have shown that such peaks 
are expected to be highly modulated, both in amplitude as well as in shape, which allows 
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Figure 12: Under the assumption of a peak at 41keV on a tungsten target, we show the regions 
accessible to xenon targets for the same dark matter stream. In regions where xenon detectors are 
sensitive, we show the position of the low energy peak in the spectrum. 

them to be distinguished from other possible origins of peaked signals, such as momentum 
dependent scattering or down-scattering of light states. 

We have applied these considerations to the CRESST commissioning run, which had 
three events clustered near 41 keV. We find that we can naturally interpret these events as 
arising from inelastic dark matter in the presence of a dark matter halo velocity structure. 
The particle physics parameters needed to explain this signal are broad and generally 
consistent with those needed to explain DAMA, although it appears the DAMA modulation 
would not arise from the same velocity structure. In some ranges of parameters, upcoming 
xenon experiments will be capable of seeing an associated peak, although for much of the 
parameter space, only tungsten targets are sensitive. If this peak is truly arising from a 
scenario such as we have described, the modulation should be easily visible in an annual 
cycle of CRESST data. 

While the situation we have considered is speculative, it stresses the need to be aware 
of such possibilities for the expected signals. The dramatic differences of such signals to 
those expected from standard elastic WIMPs in an isothermal halo can inform approaches 
to future searches. Should this scenario turn out to be realized in nature, this would allow 
for remarkable opportunities to acquire a wealth of information for both particle physics 
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and astrophysics. 
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